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We have identified a novel mechanism for the birth of Strange Nonchaotic Attractor (SNA) in a 
quasiperiodically forced Chua's circuit. In this study the amplitude of one of the external driving 
forces is considered as the control parameter. By varying this control parameter, we find that 
bubbles appear in the strands of the torus. These bubbles start to double in number as the control 
parameter is increased. On increasing the parameter continuously, successive doubling of the bubbles 
occurs, leading to the birth of SNAs. We call this mechanism as the bubble doubling mechanism. 
The formation of SNA through this bubble doubling route is confirmed numerically, using Poincare 
maps, maximal Lyapunov exponent and its variance and the distribution of finite-time Lyapunov 
exponents. Also a quantitative confirmation of the strange nonchaotic dynamics is carried out with 
the help of singular continuous spectrum analysis. 

PACS numbers: 05.45.Ac, 05.45.Pq 



I. INTRODUCTION 

Strange Nonchaotic Attractors (SNAs) are known 
to appear in quasiperiodically forced dynamical sys- 
tems. They posseses a complicated geometrical struc- 
ture, which shows them to be fractal in nature, but do not 
have any sensitive dependence on initial conditions, as 
seen from a negative maximal Lyapunov exponent. Fol- 
lowing the pioneering work by Grebogi et al. pQ, many 
researchers have found the existence of the SNAs, and 
have even classified the mechanisms for their occurrence 
in quasiperiodically forced continuous dynamical systems 
and maps. In particular the SNAs have been reported to 
arise in many physically relevant situations such as the 
quasiperiodically forced pendulum [5] , the quantum par- 
ticles in quasiperiodic potentials [3j , biological oscillators 
[4], the quasiperiodically driven Duffing- type oscillators 
[SH2], velocity dependent oscillators [10], electronic cir- 
cuits [ITHT5] , and in certain maps [TtrfPoI . Also, 
these exotic attractors were confirmed by an experiment 
consisting of a quasiperiodically forced, buckled, magne- 
toelastic ribbon |26] . in analog simulations of a multi- 
stable potential [27J , and in a neon glow discharge exper- 
iment 28 . The SNAs are also related to the Anderson 
localization in the Schrodinger equation with a quasiperi- 
odic potential 29 , 30 and they may have a practical ap- 
plication in secure communication [31-33 . 
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Broadly the mechanisms for the birth of SNAs dif- 
fer with the routes to SNA followed by the dynamical 
systems. A list of the various mechanisms pertaining 
to different routes to chaos are listed in Table - I [34] , 
Very recently, a new route namely the bubble route to 
SNAs have been reported by Senthilkumar et al. [33J 
in a quasiperiodically forced negative conductance series 
LCR circuit with a diode, wherein, one of the driving 
force is a non-sinusoidal (square wave). In this route 
bubbles appear in the strands of the torus as a function 
of the control parameter. These bubbles then grow in size 
as the control parameter is increased. Subsequently the 
strands of the bubbles start to wrinkle, giving birth to 
SNAs. The reason for this is that the quasiperiodic orbits 
become increasingly unstable in the transverse direction 
with the increase of the control parameter. This insta- 
bility is induced by the square wave type quasiperiodic 
force resulting in an increase in the size of the doubled 
strands (bubbles), followed by an enhanced wrinkling of 
these. 

In this present work, we report successively doubling of 
the bubbles route to SNAs in a quasiperiodically forced 
Chua's circuit. Here, the bubbles appear in the strands 
of the torus as a function of the control parameter, then 
grow in size and get doubled in number. After many 
successive doubling of the bubbles, SNAs are born. All 
through these, the remaining part of the strands of the 
torus which do not contain the bubbles, remain unaf- 
fected even though varying the control parameter. In 
this route no wrinkling of the bubbles are found to oc- 
cur |35j . In order to confirm this, we present here a 
detailed numerical analysis using Poincare map, Fourier 
spectrum, maximal Lyapunov exponent, its variance, fi- 
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TABLE I: Routes and mechanisms for the formation of SNAs 



Type of route 


Mechanism 


Heagy-Hammel [25] 


Collision of period-doubled 
torus with its unstable 
parent 


Gradual Fractilization [2D] 


Increased wrinkling of torus 
without any interaction with 
nearby periodic orbits 


On-off intermittency [7] 


Loss of transverse stability 
of torus 


Type-I intermittency [TJ] 


Due to saddle-node bifurca- 
tion, a torus is replaced by 
SNA 


Type-Ill intermittency H31 


Subharmonic instability 


Homoclinic collision [30] 


Homoclinic collisions of the 
quasiperiodic orbits 




Fi(t) 



F 2 (t) 



FIG. 1: ircuit realization of the quasiperiodically forced 
Chua's circuit. Here Nr is the Chua's diode. The param- 
eter values of the circuit elements are fixed as Ci = 16.8 nF, 
C 2 = 98.0 nF, L = 18.75 mH, R = 1565 fl , G a = -0.758 
mS, G b = -0.409 mS, v x = 2849.0 Hz, v 2 = 3560.0 Hz and 
F2 — 5 mV are chosen respectively. The forcing amplitude 
Fi is taken as control parameter which being varied in our 
analysis. 



nite time Lyapunov exponents and separation of near by 
trajectories in phase space. Further the singular continu- 
ous spectrum analysis is used to find scaling the exponent 
and fractal paths. 

The paper is organized as follows. In Sec. II, we 
discuss the circuit realization of the quasiperiodically 
forced Chua's circuit. In Sec. Ill, the observations of the 
bubble doubling route to SNAs are given with the help 
of Poincare cross section. In Sec. IV, we substantiate 
the observations by the maximal Lyapunov exponent 
spectrum and its variance. In Sec. V, we find the scaling 
exponent and fractal path by using singular continuous 
spectrum analysis. In Sec. VI, we discuss another 
qualitative measure of this route by separation of near 
by trajectories. In Sec. VII, the finite time Lyapunov 
exponents are used to confirm the SNA and show the 
mechanism of the bubble doubling route. Finally, in Sec. 
VIII, we summarize our results. 



II. CIRCUIT REALIZATION 

The quasiperiodically forced Chua's circuit [12] is 
shown in Fig. (1). Using Kirchoff's laws, the equations 
of this circuit derived are 



c 2 



dv, 



dt 
dvg 2 
dt 
di L 



L 



dt 



G(v C2 - v Cl ) -g(v Cl ), (la) 
G(v Cl - v C ' 2 ) + h, (lb) 
-v C2 + F\ sin(wit) + F 2 sin^i). (lc) 



Here, 



g(v Cl ) =G b v Cl +0.5(G a ~G b )[\v Cl +B p \ - \ VCl -B p \](ld) 



where vc t , vc 2 are the voltages across the capacitors C\ 
and C 2 , and is the current through inductor L respec- 
tively. Fi sin(wit) and F% sva(ui 2 t) are the external ex- 
citations having incommensurate frequencies, g{vc t ) is 
mathematical expression of Chua's diode, where G b is 
the outer slope and G a is the inner slope of the nonlinear 
curve. Bp is the break point voltage. The parameter 
values of the circuit elements are fixed as C\ — 16.8 nF, 
C 2 = 98.0 nF, L = 18.75 mH, R = 1565 fl, G a = -0.758 
mS and G b = -0.409 mS. The value of the v x {= bj x /2%) 
and 1/2 (= W2/27r) are chosen as 2849.0 Hz and 3560.0 Hz 
respectively. The forcing amplitude F 2 is fixed as 5 mV 
and the forcing amplitude F\ is assumed as the control 
parameter. In order to study the dynamics of the circuit, 
the circuit Eqs. (1) can be normalized using the follow- 
ing rescale variables and parameters. vc x — xB p , vc 2 — 
yB p , i L = B p Gz, t = C 2 t/G, G = 1/R, a = RG a , b 
= RG b , a = Ct/Cx, P = C 2 R 2 /L, fx = Fx/3/B p , f 2 = 
F 2 p/B p , ujx = GVLx/C 2 and uj 2 = Gil 2 /C 2 . 



dx 

~dt 
dy 

dt 
dz 

dt 



a(y — x — h[x)) , (2a) 
x — y + z, (2b) 
-£y + / lS m(ftit) + / 2 siii(n 2 i), (2c) 



where, 



h(x) = bx + 0.5(a - b)[\x + 1| - \x - 1| 



(2d) 



Eqs. (2) are numerically integrated using Runge - Kutta 
fourth order algorithm with step size 0.0001. The val- 
ues of the normalized parameters are a = —1.1862, b = 
-0.6400, a = 5.8333, /3 = 12.8012, fix = 2.7454, il 2 = 
3.4306 and f 2 = 0.1280. The normalized parameter fx of 
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external quasiperiodic forcing is taken as control param- iA 
eter. 



(b> 



III. BUBBLE DOUBLING ROUTE TO SNA 



To observe the emergence of SNA the amplitude of 
the external forcing fx is varied in the range fx £ 
(0.7296,0.8500). The Poincare surface of section plot of 
the two strands corresponding to period— 2 torus for the 
value of fx = 0.7296 is shown in Fig. 2(a) and the cor- 
responding power spectrum are depicted in Fig. 3a(ii). 
As the value of the parameter fx is increased further, 
bubbles appear in the strands of the torus, for a crit- 
ical value fx = 0.7744 these are shown in Fig. 2(b). 
Further, increasing the fx value from 0.7744 the bubbles 
get doubled as shown in Figs. 2(c & d). Beyond the 
value of fx = 0.7803, the bubbles get successive doubling 
Fig. 2(e) and cause the emergence of SNA. The SNA 
start to appears above the value of fx — 0.8123. It is to 
be noted that the remaining part of the strands of the 
torus remain unaffected. The Poincare surface section 
plot of SNA for the value of fx = 0.8200 is shown in Fig. 
2(f) and the corresponding power spectrum is shown in 
Fig. 3b(ii). Beyond the value of f x = 0.8200, the bi- 
furcated strands corresponding to period— 2 torus merge 
into a single strand and then transit to chaotic attractor, 
which are shown in Fig. 2(g) and Fig. 2(h) respectively. 
The power spectrum of the chaotic attractors is shown 
in Fig. 3c(ii). These observations show clearly that the 
bubbles that appear in the strands of the torus undergo 
successive doubling and cause the emergence of strange 
nonchaotic dynamics. The mechanism for this route is 
that the quasiperiodic orbit becomes unstable as a func- 
tion of control parameter and resulting the appearance 
of bubbles in certain parts of the main strands. Further 
increasing the control parameter these bubbles losses it's 
stability and gets successive doubling and leads to birth 
of SNAs. These observations point out that this route 
is significantly different from other routes known in dy- 
namical systems |34| . 
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FIG. 2: Projection of the numerically simulated Poincare sur- 
face section plots of the attractors of Eqs. (2) in the (cj> — x) 
plane with cj> modulo 2ir as a function of amplitude of the forc- 
ing fx including the transition from quasiperiodic attractors 
to SNA through bubble doubling route: (a) period— 2 torus 
for fx = 0.7296, (b) bubbled strands of period— 2 torus fx = 
0.7744, (c & d) doubled bubble strands of the period— 2 torus 
for fx = 0.7770 and fx = 0.7803, (e) bifurcated bubble, (f) 
SNA for fx — 0.8200, (g) merged two torus and (h) chaotic 
attractor for fx = 0.9475. 



IV. LYAPUNOV EXPONENT AND ITS 
VARIANCE 



1 M 



(3a) 
(3b) 



To characterize the quasiperiodic, strange nonchaotic 
and chaotic dynamics, the largest Lyapunov exponent 
and its variance are calculated as functions of fx for the 
range of 0.8000 to 0.8500. The variance(cr) of the largest 
asymptotic Lyapunov exponent (A) from finite time Lya- 
punov exponents (A) of length M, and its variance defined 



where K is the total number of finite LE (A) and A, is 
the instantaneous Lyapunov exponent. The plot of the 
largest Lyapunov exponent and its variance are shown in 
Figs. 4(a & b). It clearly demarcate the range of the 
three dynamical behaviours such as quasiperiodic, SNA 
and chaos. Below the value of fx — 0.8123 the system 
exhibits quasiperiodic behavior. In the range fx — 0.8123 
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FIG. 3: Projection of the numerically simulated attractors of 
Eqs. (2) with (i) Poincare plot with <f> modulo 2ir in the (cj>—x) 
plane and (ii) corresponding power spectrum for various value 
of fi indicating the transition from a quasiperiodic attractor 
to an SNA and chaos, (a) period— 2 torus for /i = 0.7296, 
(b) SNA for fx = 0.8200 and (c) chaos for /i = 0.9475. 



to /i = 0.8292 the largest Lyapunov exponent is negative 
and the corresponding variance is higher in this range, is 
indicate the existence of SNA. The strange nonchaotic 
dynamics becomes chaotic dynamics beyond the value of 
fi =0.8292 (after the largest Lyapunov exponent crosses 
zero axis). 
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FIG. 4: Transition from two torus to chaos through SNA. 
(a) the largest Lyapunov exponent A and (b) variance a as a 
function of parameter /i. 
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V. SINGULAR CONTINUOUS SPECTRUM 
ANALYSIS 

The singular continuous spectrum analysis, which is 
used to quantitatively confirm the strange nonchaotic na- 
ture of the dynamics, was first proposed in the investi- 
gation of model of quasiperiodic lattices and quasiperi- 
odically forced quantum systems [37] ■ In general, power 
spectra of dissipative dynamical system can be either dis- 
crete, or continuous, or a combination of both. Discrete 
spectrum are usually generated by regular motions such 
as periodic or quasiperiodic motions, where as continuous 
spectra correspond to irregular motions such as chaotic 
or random motion. A singular continuous spectrum is 
a mixture of both discrete and continuous spectra [7]. 
Using this property, we can identify whether the bubble 
doubling indeed leads to strange nonchaotic dynamics or 
not. To confirm this we compute the partial Fourier sum 



where a is proportional to the irrational driving fre- 
quency f2i and {xk} is the time series of the variable 
x of length N for the value of f x = 0.8200. When N is 
regarded as time, \X(a, N)\ 2 grows with time N as [35] 

|AT(a,A0| 2 - N» (5) 

where fj, is the scaling exponent. The time evolution 
of X(a,N) can be represented by an orbit or a walker 
in the complex plane (Re[X(a,N)],Im[X(a,N)]) and 
for a singular continuous spectrum (when the dynam- 
ics is strange) it implies that the walk on the plane 
(Re[X(a,N)],Im[X(a,N)]) will be a fractal [H>J. Fig. 
5(a) shows, the plot of N vs |A"(a,A^)| 2 for a = Qi/4, 
which has the scaling exponent fj, f=s 1.345 and the walk in 
Fig. 5(b) appears to be fractal. These results, 1 < fi < 2 
and fractal walk, strongly suggest that the dynamics is 
indeed strange and nonchaotic. 
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VI. SEPARATION OF NEAR BY POINTS 

SNAs exhibit complicated geometrical structure like 
chaotic attractors. One way to distinguish SNAs from 
chaotic attractors is to look for the sensitive dependence 
on initial conditions |10j . In order to verify this sensitive 
dependence, we analyse the separation A between two 
orbits starting from two near by initial conditions. For 
this two near by points on the attractor (%i,yi,Zi) and 
(xj,yj, Zj) are chosen and their separation 

A = V(*i-*i) 2 + (l/i-w) 2 + (*-*i) 2 (6) 



is monitored at each forward step. For regular be- 
haviour A will decay to zero as t— > oo, but for chaotic 
behaviour, A becomes irregular. In Figs. (6) a plots of 
t versus A is shown. In Fig. 6(a) shows A for the value 
of /i = 0.8200 diminishes to zero in short interval of t 
and Fig. 6(b) shows A for fx — 0.9475 sustains the irreg- 
ular variation for long time (t— > oo). Hence the former 
case clearly supports the loss of sensitive dependence on 
initial condition, which corresponds to the strange non- 
chaotic attractor, while the later one corresponds to a 
chaotic attractor. 



FIG. 6: Separation of two nearby points vs time (t). (a) 
strange nonchaotic attractor for the value of /i = 0.8200 and 
(b) chaotic attractor for the value of /i = 0.9475. 



VII. DISTRIBUTION OF FINITE TIME 
LYAPUNOV EXPONENT 

We have discussed about the bubble doubling dynam- 
ics qualitatively in section III, by using Poincare surface 
of section plots in the (<j> — x) plane, also quantitatively 
by using maximal Lyapunov exponent and sigular contin- 
uous spectrum analysis in section IV. In addition to the 
above it is also possible to distinguish torus and SNA us- 
ing distribution of finite time Lyapunov exponent |36j . It 
has been shown [14] that a trajectory on a SNA actually 
possesses positive Lyapunov exponent in finite time inter- 
vals, although the asymptotic exponent is negative. As 
a consequence, it is possible to observe different charac- 
teristic of the SNA created through different mechanisms 
[3"5] by studying the differences in the probability (P) of 
distribution of finite time Lyapunov exponent P(M, A) 
for positive and negative values. We have calculated the 
distribution of finite time Lyapunov exponent P(2000, A) 
for the attractor shown in Fig. 3b(i) and have plotted the 
same in Fig. (7). Here we find that the distribution ex- 
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FIG. 7: Distribution of finite time Lyapunov exponent is 
taken for M = 2000 for the value of /i = 0.8200 

hibits an elongated tail for the negative values, thereby 
confirming the existence of bubble doubling transition to 
SNA. This may be explained as due to the fact that in 
the bubble doubling transition, the unaffected regions of 



the strands of the period— 2 torus remain so even after 
the birth of SNA. 



VIII. SUMMARY AND CONCLUSION 

In this paper, we have reported a novel machanism for 
the birth of strange nonchaotic attractor in quasiperiod- 
ically forced Chua's circuit. We term this route as bubble 
doubling route to SNA. At first, we have presented the ap- 
pearence of bubble through Poincare surface of cross sec- 
tion and the bifurcation of the bubble by varing the con- 
trol parameter f\ in the suitable range. Further, we used 
power spectrum to qualitatively distinguish quasiperiod- 
icity, SNA and choatic attractor. Followed by this, we 
have computed Lyapunov exponent and its variance as a 
function of control parameter fx and these clearly distin- 
guish the dynamical region of quasiperiodicity, SNA and 
chaos. To confirm the presence of SNA we have com- 
puted singular continuous spectrum, which clearly shows 
the fractal dimension of the SNA and its fractal path in 
the complex plane. The separation of near by point test 
and distribution of the finite time Lyapunov exponent 
also explains the presence of SNAs. It is to be noted 
that the formation of SNA through this bubble doubling 
mechanism is significantly different from the well known 
mechanisms in the literature (Table - I). 
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